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NILPOTENT INTEGRABILITY, REDUCTION OF 
DYNAMICAL SYSTEMS AND A THIRD-ORDER 
CALOGERO-MOSER SYSTEM 

A. IBORT 1,2 , G. MARMO 3 , M. A. RODRIGUEZ 4 , AND P. TEMPESTA 2 ’ 4 


Abstract. We present an algebraic formulation of the notion of inte- 
grability of dynamical systems, based on a nilpotency property of its 
flow: it can be explicitly described as a polynomial on its evolution pa¬ 
rameter. Such a property is established in a purely geometric-algebraic 
language, in terms both of the algebra of all higher-order constants of 
the motion (named the nilpotent algebra of the dynamics), and of a 
maximal Abelian algebra of symmetries (called a Cartan subalgebra of 
the dynamics). It is shown that this notion of integrability amounts 
to the annihilator of the nilpotent algebra being contained in a Cartan 
subalgebra of the dynamics. Systems exhibiting this property will be 
said to be nilpotent integrable. 

Our notion of nilpotent integrability offers a new insight into the in¬ 
trinsic dynamical properties of a system, which is independent of any 
auxiliary geometric structure defined on its phase space. At the same 
time, it extends in a natural way the classical concept of complete inte¬ 
grability for Hamiltonian systems. 

An algebraic reduction procedure valid for nilpotent integrable sys¬ 
tems, generalizing the well-known reduction procedures for symplectic 
and/or Poisson systems on appropriate quotient spaces, is also discussed. 

In particular, it is shown that a large class of nilpotent integrable sys¬ 
tems can be obtained by reduction of higher-order free systems. The case 
of the third-order free system is analyzed and a nontrivial set of third- 
order Calogero-Moser-like nilpotent integrable equations is obtained. 
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1. Introduction 

In this paper, a new approach to the notion of integrability of dynamical 
systems that extends the standard approach a la Liouville, and also includes 
such obvious “integrable” systems like uniformly accelerated ones, will be 
discussed. For reasons that will become evident later on, this new notion of 
integrability will be called nilpotent integrability. 

One of the main features of nilpotent integrability is that it does not 
require any auxiliary underlying geometrical structure for its formulation. 
In this perspective, nilpotent integrability becomes an intrinsic property of 
a dynamical system, and the specific geometry associated to it comes later. 
Moreover, we shall show that nilpotent integrability is consistent with the 
different procedures of reduction available for dynamical systems. They will 
allow to create new nilpotent integrable systems out of simple examples. 

The main motivation behind the idea that ‘integrability’, as a property of 
a dynamical system, must be prior to any geometry associated to it, comes 
from a most pragmatic approach to integrability. For instance, if we demand 
that a system is ‘integrable’ if its flow could be explicitly exhibited by using 
a finite number of quadratures, there is no obvious geometrical content in 
such a definition. 

As a logical consequence of this observation, it is natural to try to refor¬ 
mulate the well established notion of complete integrability of Hamiltonian 
systems without appealing to any auxiliary geometry, such as a symplectic or 
a Poisson structure (see, for instance, |Mm85| and references therein). The 
consistency of these ideas with general reduction procedures would reflect 
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the fact that quadratures and elimination of variables should be compatible 
with each other. 

To clarify the scope of the ideas leading to the notion of nilpotent integra- 
bility, let us consider first the instance of a complete integrable Hamiltonian 
system system T from the standard perspective. Thus, if T represents a com¬ 
plete integrable Hamiltonian system on a symplectic manifold, introducing 
action-angle variables hj, , j = 1 ,... ,n, the vector field describing T will 
look like 


Q 

v ’ dOi 

The explicit expression of its flow ipt, given by 


hj o (p = hj , o ipt = + v-* (h)t , 


shows that the flow ipt, described in action-angle variables, is affine in t. 
More important, it is trivial to check that the derivatives along the dynamics 
of the chosen coordinates are given by 

(1) L r e j = u j (h)- C v hj = 0. 


The previous equations can also be interpreted by saying that the repre¬ 
sentation in action-angle variables of the derivation T is nilpotent of order 
2. It is relevant to point out that Eqns.Q can be stated as a definition of 
complete integrability without any reference to a symplectic structure. 

Let us consider now simple examples of systems which are ‘integrable’ 
in an obvious sense but that do not fit naturally in the previous scheme. 
For instance, let us consider a uniformly accelerated system in M 3 : r = 0, 
r € M. 

The equations of motion of such a system are described by the vector field 
on (r, v, a) 6 T 2 M 3 


d d 

r = v- 7T +a ' o - ’ 

OT OV 

with a a constant vector. Its flow tpt has the explicit polynomial expression 


r o <p t = r + vt -I- -at 


Notice that the acceleration a is obviously a constant of the motion, while v 
and r are higher-order constants of the motion, that isv = r, a = v = r and 
a = 0. Similarly, the kinetic energy E\ = ^v 2 is not anymore a constant of 
the motion. However, 


E- 2 . — Ei — v • a, — E 2 — a 2 , F 3 — 0 , 


i.e., E 3 is a constant of the motion, while E 2 , E\ are higher-order constants 
of the motion (of orders 2 and 3 respectively). The same can be argued 
with respect to the tower of (higher-order) constants of the motion derived 
from the standard angular momentum: Li = r A v, L 2 = Li = r A a, 
L 3 = L 2 = v A a, L 3 = 0. 
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The way how the previous ideas fit with the notion of reduction can be 
nicely described by considering the basic example of the completely inte¬ 
grate free system r = 0 in R 3 . It is easy to construct other, perhaps more 
interesting, completely integrate systems by reducing it using its integrals 
of motion. For instance, by using the standard polar decomposition of r, 

r = rn, n • n = 1 ; 

r = rn + rri , n • ri = 0 ; 

r = rn + 2 rri + rri, n • ri = — ri 2 , 

we may write its constants of the motion E\ and Li as 

Ei = - (r + r n ) , l 1 =L, 1 =rn. 

Then, once we restrict to the submanifold 

£(fc, a) = {(r, v) € TR 3 | al\ + 2(1 — a)E\ = k, 0 < a < 1} , 

since for our system we have r = rri 2 , we get the following family of 
completely (and explicitly) integrate reduced second-order systems in one- 
dimension: 

al 2 + 2(1 — a) Ei — (1 — a)r 2 
r(ar 2 + (1 — a)) 

It is immediate to generalize the previous construction to the system r = 0. 
Now we have the additional relations 

r = r n + 3rri + 3rri + rri , n • ri = —3ri • ri . 

Moreover, we have r = 3rri 2 — rn • ri. The towers of higher-order constants 
of the motion Ei, E- 2 , E 3 and Li, L 2 , L 3 before dehned, allow us to write e.g. 
E ‘2 = ri • ri = rr + rrri 2 + r 2 ri • ri. Hence we get the explicitly integrable third 
order system: 

... = 3 E 2 - rr = a 2 t + c. 

r 

In other words, the reduction of systems like r = 0, possessing polynomial 
flows, gives rise to systems that can be explicitly integrated and possess 
polynomial flows as well. Additional examples of this kind will be discussed 
later on. 

Therefore, the main idea behind the notion of nilpotent integrability is to 
take advantage not only of the constants of the motion, but also of the whole 
family of higher-order integrals of motion associated with the system. If this 
family (that will be called the nilpotent algebra of the dynamical system) 
is large enough, then the flow of the system will be shown to be explicitly 
described, in terms of the integrals, as a polynomial function of time. In such 
a case, the system will be said to be strictly nilpotent integrable. Systems 
obtained by reduction of the system r = 0 will provide examples of the 
strictly nilpotent case. 

However, it could also happen that the nilpotent algebra associated with 
a given system would not be large enough, although there exists an Abelian 
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algebra of symmetries of the dynamics large enough to complement it. For 
instance, in the case of completely integrable systems such Abelian algebra 
of symmetries is the algebra of Hamiltonian vector fields defined by the ac¬ 
tion variables, whose flow parameters are the “angle” variables 0 3 . Such 
systems will be called the completely nilpotent integrable ones. A partic¬ 
ular instance of them are the standard completely integrable Hamiltonian 
systems to which we have referred to. 

Due to the paramount role played in this analysis by the algebraic struc¬ 
ture of the families of higher-order constants of the motion and various Lie 
algebras of symmetries of the system, we are led to develop the theory in 
an algebraic framework, finding also inspiration in the notion of algebraic 
integrability introduced in I,a93 . 

It is important to mention that the algebraic framework also allows to 
extend in a meaningful way the notion of ‘integrability” to quantum sys¬ 
tems (see for instance EH], fFaTTll . and the discussion in fT£H| . Chaps. 
7 and 8). In fact, the algebraic approach to the description of physical sys¬ 
tems has gained a significant weight since the discovery of quantum group 
symmetries in physical systems. More recent developments in quantum grav¬ 
ity (see for instance |La97| . [ VaOBj . Il ia 10 1 point out the need for a non- 
commutative description of space-time, hence emphasizing the algebraic ap¬ 
proach to their study. In all these contexts, the algebraic view-point instead 
of a set-theoretical description of physical systems is mandatory. 

Hence, it is all but natural to pursue a foundation of the theory of dy¬ 
namical systems from this point of view. Even more, in doing so, as the 
following discussion will show, new insight and ideas can emerge naturally. 

Consequently, we shall introduce the notions of strict and complete nilpo¬ 
tent integrability of a dynamical system T by identifying it with a derivation 
of an algebra (of an appropriate class) T. We will prove that, under adequate 
conditions, the flow of such system can be represented as a polynomial func¬ 
tion of time. As was mentioned above, our theory applies to a large class 
of dynamical systems, possessing no obvious geometrical structure, like a 
symplectic, or Poisson one, etc. The algebraic framework will also make 
transparent the consistency of generalized reduction procedures (i.e. reduc¬ 
tion with respect to subalgebras and ideals compatible with the dynamics) 
with the notion of nilpotent integrability. Hence we shall prove that, under 
appropriate hypotheses, if a system is nilpotent integrable its reductions are 
nilpotent integrable too. 

The formalism of differentiable spaces is the language we have chosen 
to deal with the algebraic formulation of dynamics. Differentiable spaces 
are one of the most useful extensions of the notion of differentiable mani¬ 
folds that permit a simple algebraic description. In fact, local models for 
differentiable spaces are just differentiable algebras, i.e., quotients of alge¬ 
bras of smooth functions on open sets in finite dimensional Euclidean spaces 
equipped with the strong Whitney topology. The use of differentiable spaces 
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will also allow to keep the formalism close to the standard description of dy¬ 
namical systems on manifolds. Another advantage of this approach is that 
a large class of mild singularities are automatically dealt with within this 
formalism. Also, some technical difficulties in dealing with local charts in 
the standard formalism of differentiable manifolds are overlooked. It is also 
useful to mention here that the present algebraic approach in principle can 
be extended to the case of dynamical systems in the context of Galois dif¬ 
ferential algebras along the lines discussed for instance in hee]. 

The structure of the paper is the following. The first subsections of Sec¬ 
tion [2] will be devoted to review the standard notions of dynamical systems 
both in the language of smooth manifolds and of differentiable spaces, mak¬ 
ing the transition as smooth as possible. In Section [3j we shall discuss the 
fundamental structures attached to a dynamical system, i.e., its algebra of 
symmetries and its nilpotent algebra, and we will introduce the fundamental 
notion of nilpotent integrability. In Section [4j a theory of generalized reduc¬ 
tions will be formulated in the framework of nilpotent integrability. Finally, 
in Section [5] a family of examples leading to a Calogero-Moser higher-order 
nilpotent integrable system will be discussed. 


2. Algebras, derivations and dynamical systems 


2.1. Dynamical systems in smooth manifolds. Let M be a smooth 
paracompact differentiable manifold of dimension n. Local coordinates on 
M will be denoted by x l , i = 1 ,...,n. Let F(M) be the commutative 
algebra of smooth functions on M. Hereafter we will restrict to real valued 
functions (although in some applications complex valued functions may also 
be relevant). 

Dynamics on M will be described by a vector field T, which is a derivation 
of the algebra F(M), i.e. a linear map on F(M) such that T(fg) = T(f)g + 
fT(g). The space of derivations of the algebra F(M), denoted in what 
follows by Der(J r (Ad)), is a module over the algebra F(M). Moreover, it 
converts into a Lie algebra if equipped with the Lie bracket 

[Y 1 ,Y 2 ](f) = Y 1 (Y 2 (f))-Y 2 (Y 1 (f)), VYi,y 2 €Der (F(M)), f 6 F(M). 

The dynamics represented by T will be described in two possible ways, 
a) As a system of autonomous first order differential equations, that in local 
coordinates x l takes the form: 


( 2 ) 


dF_ 

dt 


I», i = l,...,m, 


where the vector field locally reads T = T l (x)d/dx 1 . 
b) Equivalently, as an evolution equation on F(M ): 


dj_ 

dt 


(3) 


F(/)i / 6 F(M). 
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The local flow of T will be denoted by pt and satisfies: 


(4) 


d<pt 

dt 


r o p t . 


If the vector field is complete, the flow pt defines a one-parameter group 
of automorphisms of the algebra J~(M) (a one-parameter family of auto¬ 
morphisms for non-autonomous systems). Alternatively, we can integrate 
formally eq. ©> to obtain: 

°° j-Tl 

(5) P* t f = f°Pt = exp(tr)/ = ^— T n (f). 


2.2. Algebras of functions and differentiable algebras. The algebra 
J-(M) of smooth functions on the manifold M is equipped with the topology 
of uniform convergence on compact sets of the function and their deriva¬ 
tives to order r ( r = 1,2,... ,00). However, for the purpose of extending 
the standard notions of differential calculus, it is necessary to consider the 
strong or Whitney topology on IF(M), that is the topology defined by the 
countable family of seminorms Pik(f) = where {Kf\ is 

a subcovering by compact sets of a small enough open covering \Uf\ of M. 
Thus J~(M) becomes a Frechet topological algebra, i.e., a locally multiplica¬ 
tive convex algebra which is metrizable and complete, where multiplicative 
convex means that the topology is defined by a countable family of submul- 
tiplicative seminorms. 

The algebra C°°(M n ) is an instance of a Frechet topological algebra and 
the model for the notion of a differentiable algebra [Ma 66 | (see also |Na03] 
for a recent account of the theory of differentiable algebras). 

Definition 1 . A differentiable algebra is a real Frechet algebra isomorphic 
to C ,0 °(M m )/ l 7, where J is a closed ideal ofC°°(M m ). 

Notice that the class of differentiable algebras is larger than the class of 
algebras of differentiable functions on a manifold (for instance C°°(M)/(x 2 ) 
is not an algebra of differentiable functions on a smooth manifold). 

Because of Whitney’s embedding theorem, any smooth paracompact man¬ 
ifold M can be embedded into M m , hence J r (M) = C°°(W n )/J', where J 
is the closed ideal of smooth function on M m vanishing on M (as a closed 
submanifold of it). Hence J- = C°°(M ) is a differentiable algebra too. 

Definition 2. Given a real topological algebra T , we define its spectrum 
Spec R (J 7 ) as the space of continuous M.-morphisms ip: T — M equipped 
with the natural topology. 

This space can be identified with the space of real maximal ideals of T. 

When M is a Hausdorff smooth manifold satisfying the second countabil¬ 
ity axiom, then Spec R (J r (M)) = M. Given an element / e T and a point 
x E Spec^J 7 we define the GePfand transform f(x) = x(f). Thus, any / 
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defines a function on Spec R J-\ If we consider then / coincides with 

/• 

The spectrum of a differentiable algebra is the local model for a class of 
spaces called C°°-spaces or differentiable spaces for short, that extend in a 
natural way the notion of smooth manifolds (see for instance [Na03j and 
references therein for a modern account of the foundations of the theory of 
differentiable spaces). 

Differentiable spaces are ringed spaces over a sheaf of differentiable alge¬ 
bras. They allow to work with a large class of topological spaces, which share 
with smooth manifolds most of their theoretical formulation and the differ¬ 
ential calculus attached to them. At the same time, differentiable spaces 
can exhibit singularities and other structures typically arising in the stan¬ 
dard processes of reduction of dynamical systems. As was commented in 
the introduction, this is the reason why we have chosen such a structure 
to introduce the notion of nilpotent integrability. Actually, we will use in 
this paper the “local” formulation of the theory only, i.e., the notion of dif¬ 
ferentiable algebras. However, once formulated in the case of differentiable 
algebras, nilpotent integrability can be extended naturally to the category 
of differentiable spaces. 

Differentiable spaces defined as the real spectrum of a differentiable al¬ 
gebra are called affine differentiable spaces. It can be proved that a dif¬ 
ferentiable space is affine if and only if it is separated. This means that at 
each point there is a finite local family of functions separating infinitely near 
points (see later, Def. 2), with bounded dimension, and its topology has a 
countable basis (see the “Embedding Theorem” at |Na03j . p. 67). Thus, 
any compact separated differentiable space is affine and a standard smooth 
manifold can be characterized as a separated, finite-dimensional differen¬ 
tiable space with a countable basis for its topology. 

Given a differentiable algebra J 7 and a point x 6 Spec^J 7 , let be the 
maximal ideal of all / £ J such that f{x) = 0. Also, let T x := T mx be 
the localization of J 7 with respect to the multiplicative system S x defined 
by rriz, namely the set of all elements f ^ m x . 

Definition 3. We shall call the elements on T x the germs of elements f of 
J 7 at x. 

They can be considered as residue classes of elements on T with respect 
to the ideal of elements of J 7 vanishing on open neighborhoods of x\ we shall 
denote them by [f]x- 

For a differentiable algebra J 7 , the space of its derivations Der(J r ) can 
be equipped with its canonical Lie algebra structure [■, •]. Such a space is 
a ^-module. The space Der(J 7 )* of J-"-valued J r -linear maps on Der(J r ) is 
the space of 1-forms of the differentiable algebra J 7 . The space Der(J 7 )* is 
again a J r -module; we will denote it by D 1 (J 7 ) (or simply by D 1 for short). 

Let U be an open set in a differentiable space X. All objects considered 
so far can be localized to U in a standard way. Therefore, we can consider 
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the differentiable algebras Xjj, DepJ 7 );/, etc. Moreover, differentiable 
spaces admit partitions of unity (see [ Na03l p. 52). 

The differential map d: T —>■ Q 1 is defined as df(Y ) = Y(f) for any 
Y G Derp 7 ). In the case of T = the space fi 1 is just the module 

of smooth 1-forms II 1 (M) on the manifold M. 

Alternatively, we can associate with any point x G Spec^J 7 the linear 
space m x /m 2 x denoted by T*X where X = SpecjjJ 7 . Again, the space T*X is 
the localization of the module II 1 with respect to the multiplicative system 
defined by m x . Then given / G m x , we have that df(x) := [df] x can be 
identified with the class of / in m x /m x . In the same vein it is possible to 
show that the dual space to T*X, namely the tangent space T x X to the 
differentiable space X at x, is just the space of derivations of the localized 
algebra T x . 

2.3. Generating sets. 

Definition 4. Given a differentiable algebra T, we shall say that a family 
of functions g \,..., g n separate points infinitely near to x G X = SpecR.7 7 if 
the differentials dg±(x),... ,dg n (x) generate T* X as a linear space. 

Equivalently, the differentials dgk{x ) of the family separate elements v x in 
T x X , i.e., for any v x 7 ^ w x there is gu such that v x (gk) 7 ^ w x {gk)- 

Definition 5. We say that the family of functions Q = {g a } of the differ¬ 
entiable algebra T is a differential generating set if they separate points in 
X = SpecR^ 7 and if for any x G X, there is a finite subfamily {g ai ,... ,g an } 
of Q that separate infinitely near points to x. 

Notice that if Q is a differential generating set for the differentiable algebra 
J 7 , then for each pair of derivations Y\ and Y% there exists at least one element 
g a eQ such that Y^ga) f Y 2 (g a ). 

When the previous condition holds, we say that the functions g a sepa¬ 
rate derivations. It is not hard to see that a large class of manifolds, for 
instance C°° paracompact second countable manifolds, possesses generating 
sets. Also, if Q is a generating set for the algebra T as an associative algebra, 
then it is a differential generating set for T. 

Since derivations can be localized, generating sets also separate local 
derivations. This fact is used to prove the following important property 
of differential generating sets, which justifies their name. 

Lemma 1. Let M be a smooth manifold. A set of functions Q = {g a } of 
X{M) is a differential generating set iff the set of 1-forms d Q = {dg Q } gen¬ 
erates the algebraic dual of the space of derivations Derp 7 )*, which coincides 
with the space fl 1 (M) of smooth 1-forms on M. 

Proof. Consider first an open set U contained in the domain of a local chart 
of M. If there were a local 1-form crjj that could not be written as ajj = 
<ru,adg a , then the span of dg a would be a proper subspace of T*U and it 
would exist a vector field X lying in the annihilator of such subspace. Hence 
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dg a (X) = 0 for all a. This argument can be made global by using partitions 
of the unity: for any 1-form a on M it must exist a family of functions a a 


with compact support on M such that a = (J a dg a . 


□ 


In the subsequent analysis, we shall assume the existence of a differential 
generating set Q for T. Let us mention another direct consequence of the 
properties of a generating set Q. Since locally the differentials dg a generate 
T*X with X the spectrum of J~, we can extract a subset dg ai that is locally 
independent, i.e. such that dg ai A ■ ■ ■ A dg am / 0 on a open neighborhood of 
any given point. Therefore, we will say that the functions g ak define a local 
coordinate system. Later on we will use this fact to write explicit formulas 
in terms of subsets of generating sets. 

We will introduce a notion of finiteness for algebras T which is funda¬ 
mental in order to define our concept of integrability. 

Definition 6. The differentiable algebra J- is said to be of finite type if it 
admits a finite differential generating set Q = {g\, ..., gjv}, N E N. 

This more restrictive condition is satisfied for instance if the manifold is 
compact or of finite type. Indeed, in such case it can be embedded into a 
finite dimensional Euclidean space whose coordinate functions, restricted to 
the embedded manifold, provide a differential generating set. As a conse¬ 
quence of the previous discussion, it can be shown that if M is a smooth 
manifold, a differential generating set for the differentiable algebra J~(M) 
provides a set of local coordinate systems, i.e., an atlas for the manifold 
M, by restricting to small enough open sets and shieving out dependent 
functions. 

2.4. Derivations and their flows on differentiable algebras. Let T 

be a derivation of the differentiable algebra of finite type J 7 , hence T = 
C°° (M n )/ fj. Denoting by X, as before, the real spectrum of J 7 , for each 
i£l the derivation T defines an element T x E T x X. Hence, X is a closed 
differentiable subspace of M n , and the canonical injection i: X —»-M n maps 
T x to a tangent vector i*T x E Moreover, we can extend the vector 

field i*T along i(X) to a vector field T in M n . Let fit be the flow of T. By 
construction fit leaves i(X) invariant. Let us denote by tpt the restriction 
of fit to X (that always exists because of the universal property of closed 
differentiable subspaces (Na03j p. 60). We will denote by <pt the flow of the 
derivation T. The flow (ft will act on elements / E J- as 


= + f + J = f- 


Also, the flow (ft can be integrated formally by using a close analog to 
formula f[5j): 



(6) 
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We have proved that not only the derivations on differentiable algebras can 
be extended to derivations on C°°(M n ), but also that they are continuous 
maps with respect to the canonical topology on T. Indeed, they are the 
quotient of derivations on C°°(M Tl ) preserving the closed ideal J, and are 
continuous. 

Finally, we observe that formulas similar to Eqs. ©-© hold for deriva¬ 
tions on differentiable algebras: 


(7) I -IW ter, 

with the local flow tpt satisfying 


( 8 ) 


d(ft 

dt 


r o (p t . 


3. The algebraic formulation of nilpotent integrability 


3.1. Infinitesimal symmetries. 


Definition 7. Let C be a Lie algebra. Given a subset S C C, the space of 
infinitesimal symmetries of S is its commutant S', i.e. 

S , = {£g£| [f,x] = 0 ,VxgS}. 

In particular, if C = Derp 7 ) is the Lie algebra of derivations of a differen¬ 
tiable algebra T, and S = {r} ; with V e DepJ 7 ), the space of infinitesimal 
symmetries ofT is the commutant of S = {T}, i.e., the set of derivations £ 
of T such that [r,£] = 0. 

The space of infinitesimal symmetries of any subset S C £ is a Lie subal¬ 
gebra of C. 

We define recursively the nth commutant of S as: 

(9) S' = {£ G £ | [f,x] = 0,V x e S}, S( fc+1 ) = (£(*))', k> 1. 

Notice that by dehnition S C S" where S" is called the bicommutant of S. 
The following result holds. 


Lemma 2 (Stability Lemma). Let C be a Lie algebra and S C C a subset. 
Then: 

(10) ScS" = S (4) =S (6) = ••• , S' = S"' = S (5) = ••• . 

Moreover, if S is Abelian, we have S C S” C S' and the bicommutant S" is 
an Abelian Lie subalgebra of L. 


Proof. Observe that if *Si C ^2, then S' 2 C 5J. Because S C 5", we have 
S'" C S'. On the other hand, since for any set S C S" we have that 
S' C (S')", we infer that S' = S’" and relations (10) follow. 

If S is Abelian, i.e., [x, y\ = 0 for all x,y G S, then S C S' , hence S" C S'. 
Moreover if £, ( G S" C S', then ( 6 S', and [£,£]= 0, which shows that S" 
is Abelian. □ 
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Definition 8 . Let T be a derivation of the differentiable algebra T. We 
shall say that T is quasi-periodic if the commutant 5p is finite dimensional. 

Notice that if T is a vector field on a compact manifold M, then T is 
quasi-periodic if the flow of T generates a finite dimensional torus inside the 
group of diffeomorphisms of M whose Lie algebra is the commutant 5p. 

However, we will be mainly interested in the “Abelian part” of the algebra 
of symmetries of T, that is we will consider a maximal Abelian subalgebra 
of S' containing S. When the subset S is just T, we shall denote such a 
maximal Abelian subalgebra by f)r- Hence, hr is a Cartan subalgebra of 
S' containing T only; consequently, it will be called a Cartan subalgebra 
associated to the dynamics T. 

More can be said when the Lie algebra £ is represented in terms of an 
algebra of derivations of a given algebra T (or of automorphisms of an 
^-module). Let p: £—^Der(J r ) be a morphism of Lie algebras. To each el¬ 
ement £ £ £ we associate a derivation p(£), satisfying p([£, £]) = [p(£)> /°(C)]- 
We will analyse this situation in the following sections. 

3.2. Constants of the motion. 

Definition 9. Given the dynamics T on T, the subalgebra of constants of 
the motion ofT, denoted by C(T) (or simply by C) is defined by 

(11) c(r) = {/ g 7 1 r(/) = o}. 

More generally, as in the previous section we can consider a Lie algebra 
£ represented by derivations of the algebra T (in particular we can consider 
the Lie algebra Der(J r ) itself with the tautological representation). Given 
any subset S £ £, we define the subalgebra of its constants of the motion 
C(S) by 

(12) C(S) = {/ £ ^ | p{x)(f) = 0,V x £ 5}. 

We will use a notation reminiscent of the one used in the previous section by 
denoting the subalgebra of constants of the motion by S' p C T and calling it 
the “commutant” of S in T with respect to the representation p. Similarly, 
we can define the bicommutant of S in T as the Lie subalgebra of £, denoted 
as S" p . of elements £ £ £ such that />(£)(/) = 0 for all / £ C(S) = S' p C T . 

(p) 

We define recursively ST C T as follows: 

(13) = S' p = C(S), S^=S'T 

s (2k+i) = {/ £ J 7 | p(£,)(f) = 0 , £ S^} , fc > 1 

S^ = {x £ £ | p{f)(f) = 0, V/ £ 5^ 2fc_1 ) ,k >2 .} 

Then, an argument similar to that of Lemma [2] gives us another stability 
result. 

Lemma 3. The following relations hold 

(14) 5c5;=5} 4 ) = -, C{S) = S' p = s;' = •••. 
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Moreover, if S is an Abelian subset of C, then S" is the minimal Abelian 
subalgebra of C which contains S. 

Observe that the kernel of any derivation of a given algebra is a subalgebra 
of it. Since C is a subalgebra of F, if M. is a T-module, M. is a C-module 
too. In particular, F is a C— module. The space of C-linear maps from F to 
C is called the space of sections of F over C (see later for the set theoretical 
interpretation of such set). 

If F is the differentiable algebra of smooth functions on a manifold M , 
then the algebra of constants of motion for the derivation T is closed in 
F with respect to its canonical topology, because T is a continuous map. 
However, it is not true in general that a closed subalgebra of a differentiable 
algebra is a differentiable algebra. In what follows, we shall assume that 
the algebra of constants of motion for T is a differentiable algebra. This 
assumption is actually satisfied by a large family of derivations, as we will 
see in several examples. 

Let C = Spec K (C) denote the real spectrum of the differentiable algebra 
C. Hence the canonical inclusion morphism i: C —»- F induces a continuous 
map 7r: X — ^C, such that i = tt* . The injectivity of i implies that the map 
7 T is a surjective projection map. Let c G C be a point in the real spectrum 
of C. Then F/Fm c is the differentiable algebra of the fibre 7r _1 (c) of tt over 
the point c. Notice that the flow of T leaves invariant the fibres of n and 
projects to the null derivation on C. 

3.3. The nilpotent algebra of a derivation. We are able now to define 
the sets of higher-order constants of the motion or generalized constants of 
the motion. 

Definition 10. We shall say that f G F is a constant of the motion of 
order k if 

(15) ££(/) = 0 and C^f / 0. 

The number k will be called the nilpotency index of f with respect to T. The 
set of higher order constants of the motion of V up to the index k will be 
denoted by 

C k (T) = {f €F\tfKf) = 0}. 

With this notation, C(T) = C 1 (T). Observe that the sets C k (T) are closed 
linear spaces such that C fc (T) C C ; (T) for all k < l. Consequently, we can 
construct the following filtration: 


(16) C(T) =C 1 (T) C ••• CC fc (T) C ••• C 

Lemma 4. The spaces C fc (T) are C -modules. 
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Proof. By using the binomial expansion of the derivation £p acting on the 
product of two functions, we have 

(17) 4(M = e0W(sI = » 

for any / G C and g E C k . It is also clear that C l v C k C C k ~ l for all k > l and 
C l r C k = 0 for all l > k. □ 

A central definition for the present theory is the following. 


Definition 11. Given a dynamical system T in the differentiable algebra 
F, the nilpotent subspace of V is the closure of the union of all generalized 
constants of the motion. We will denote such set as Af(T): 

A/"(r) = {f eF\3ken,c k f = o}. 

The spaces C k in general are not subalgebras, contrarily to the case of AT. 
Lemma 5. The space AT is a closed subalgebra of F. 


Proof. Let /, g E Af, with nilpotency indices k\ and kp respectively. Then 
by taking p > k\ + all terms in the binomial expansion (17) of £p(/<?) 
have the form with either r or s > maxjfci, Hence, all 

terms in the expansion must vanish. □ 


Thus we have the following structure: 

IcCcAcT. 


As in the case of the algebra of constants of the motion, we shall assume 
that the subalgebra AT is a differentiable algebra. This assumption is satisfied 
by a large class of examples and applications, as will be discussed later on. 
Under this assumption, the real spectrum of Af is an affine differentiable 
space, denoted by N = Spec K (AA). Let us denote by <f>: X—*-N the canonical 
projection from X to N. There is another natural projection g: N —^ C 
that makes commutative the following diagram, i.e., such that ir = g o <f>. 


X —N 



v 

C 


Given a subset 77 C J 7 , by analogy with the definition of the commutant 
of a subset S of the Lie algebra £, we can define its commutant with respect 
to the representation p of £, i.e. we introduce VJ p = {x G C \ p(x)(f) = 
0, V/ G 77}. This notion is consistent with the definition of the bicommutant 
of S C £ given in Eq. (13). In particular, if we consider Af C F, and 
L = Der(J r ), p = id, then we put Af' = Af', with 


U' = {Y G Der(^) | Y(g) = 0, V# G Af} . 







NILPOTENT INTEGRABILITY AND REDUCTION 


15 


3.4. The flow of a regular derivation on a differentiable algebra. In 

the rest of this section, we shall assume that all of the algebras appearing 
on the discussion are of finite type. Under this assumption, the filtration 
given in Eq. (16) necessarily stabilizes, that is: 

3uEN s.t. C "- 1 ± C v = C v+1 = C v+2 = ... 


The index v will be called the index of nilpotency of T. 

We can construct a generating set for T as follows. Let {h a }, a = 
l,...,ri be a generating set for C. The C-modules A f k = C k /C k ~ l for 
k > 1 are finitely generated. For the quotients J\f k = C k /C k ~ l , k > 1 
we will choose a generating set of the form + C k ~ l (a*, = 1 ,... , rfc), 
such that T(g^) = ga k ^ ■ Notice that r(C fc ) C C fc_1 . Hence the family 
{g^\ , g $,..., g ^,..., g\ v } } is a generating set of N as a C-module. 

Definition 12. Given an element g^ E C k , the elements g = T(g ^ k ' ) ),..., 
g ( ' 1) = T(gW), h = g^ = T(g^), will be called the descendants of g^ k \ The 
set of the descendants of an element g E Af will be called the tower of higher- 
order constants of the motion defined by g and denoted by T(g). 

Similarly, given an element h E C, consider the elements g E Af such that 
T k ^ 1 (g) = h for some k. We will call them the ascendants of h. Given h, 
the set of its ascendants will be called the tower of higher-order constants of 
the motion defined by h. 


Hence the generating family constructed before decomposes as the union 
of the towers of higher-order constants of the motion defined by the contants 
of the motion h a . See Fig. [l]for a pictorial description of the structure of 

AT. 

Finally, to complete a generating set for T ’, we will extend the previ¬ 
ous family by adding a family of functions {/i,..., f n }. We will call the 
generating set so constructed a V-adapted generating set of the algebra T. 

As proved in Lemma [lj for any point x E X = Spec R (J r ), the differ¬ 
entials at x of the elements of a T-adapted generating set span the whole 
cotangent bundle TfX. Consequently we can extract from them a set of 
“local coordinates” for A in a neighborhood of x. A local coordinate set 
around the point x has the form {h\ ,..., hp , g\ 2 \ ..., g[ f'*,..., ,..., glf^ , 

1 J- 7*2 1 T u 

/i,..., f n i }. We shall call such a local coordinate system a T-adapted local 
coordinate system. 

The number of elements r^, r r 2 , etc. on each class of an adapted local 
coordinate system is locally constant, and therefore is constant on each 
connected component of X. Notice that r , 1 +r , 2 + - ■ ■-\-r , u + n' = m. = dimA, 

In what follows, we will use the same notation both for the number r\ 
of generators of C and for the number r[ of constants of the motion of an 
adapted coordinate system, and similarly for r 2 and r 2 , etc., even though on 
each class of coordinates of a T-adpated coordinate system r[ < r\, r 2 < r' 2 , 
etc. With this convention, we have that r\ + r 2 + • • • + r u + n = m. = dimAL 
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— ha j 


9a ] ° <Pt 

= 9a ] + h a t , 


9a ] ° Pt 

= 9a ] + 9 { a ] t + \h a t 2 , 

(19) 

9a ] ° <Pt 

= 9 ( a k) +9 ( a~ l) t + --- + 


kl 


3.5. Strict nilpotent integrability. We propose several definitions of nilpo- 
tent integrability. The first one is directly tied to the polynomial description 
of the flow associated with a dynamical system T, according to Eqs. (19). 


Definition 13. A dynamical system T on a differentiable algebra T is 
strictly nilpotent integrable if there exists a family of higher-order constants 
of the motion Q = {hi ,..., h r , g±,... ,g s } C M that separates points and 
infinitely near points of X = Spec^J 7 ). 


Notice that the algebra generated by Q is dense in T because it separates 
points and infinitely near points, hence in such a case A f = J-. Then the 
flow of r has an explicit description as a polynomial of the form given by 
Eqs. (19). In fact, because of the previous discussion, we can choose a local 
coordinate system around each point of the form 


{h ±,..., h ri , g ^ 


( 2 ) 


q {2) 

5 Ur2 ’ ' 


(s) 

, 9 1 , 


■ dfl 


and the flow of the dynamics T is polynomial in these coordinates. In this 
situation, we will also say that the dynamics T is nilpotent integrable of 
order s. 

Unfortunately, the above definition of nilpotent integrability, although 
it provides a direct grasp of its meaning, is hard to check. It would be 
desirable to have a definition based on the nilpotent subalgebra A f only. As 
was already pointed out, constant functions (which are always contained 
in A f) fuzz the picture. It is more appropriate to consider the algebras of 
functions vanishing at a given point xq G M. Observe that for any h a G AT 
the function h a (x ) — h Q (x o) vanishes at xq. Then we have the following 
main result. 


Theorem 1. A regular dynamics T on a differentiable algebra whose real 
spectrum is formally smooth is nilpotent integrable if and only if J-dAf = fl 1 . 

Proof. It is clear that if T is nilpotent integrable, then A f = J~, hence ob¬ 
viously FdM = Q 1 . Conversely, if FdM = fi 1 , since fi 1 defines a locally 
free sheaf over X of bounded rank (because the spectrum of T is formally 
smooth), then there exists a finite family of elements j3i,..., f3j\r generating 
it. Notice that /3 k = f k idg k i, k = 1,..., N, l = 1,... r + s, with g ki G A f. 
Consider the family Q formed by the elements g k i- Then the family of 1— 
forms dg k i generate fl 1 over T. At the same time, the localization d x g k i of 
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the 1-forms to x E X generates T*X for any x E X = Spec R (J 7 ), hence 
they separate infinitely near points. 

Moreover, let x E U with U an open set containing a given point x in the 
real spectrum of F, then the localization djjg^i of the family of 1-forms dg^i 
to the open set U generates the localization £r(U) of the sheaf of 1-forms on 
X. Consequently, we can extract a family that is independent because the 
sheaf is locally free. Let us denote by dg a such a family; then there exists 
functions g a on some open set U in M n where C °°/ J is a representation of the 
differentiable algebra F, such that U C U and g a +J = g a - The spectrum of 
an affine differentiable space is separated ( |Na03j| . p. 67). Therefore, given 
points x / y we can find open neighborhoods U and V of x and y respectively 
such that [/(1L = 0. Then the local functions g a constructed above on U 
and V provide a function g on Q separating x and y. Hence Q C AT separates 
points and infinitely near points and T is nilpotent integrable. □ 

3.6. Complete nilpotent integrability. In the analysis of integrability, 
a preeminent role is played by completely integrable systems. We have seen 
already that angle coordinates can be treated as second-order constants of 
the motion. However, in general, such coordinates do not define global 
functions, i.e., they are not defined on the differentiable algebra F. The 
explanation for this, as was already pointed out in the Introduction, is that 
angle variables should be considered in terms of an Abelian algebra of sym¬ 
metries of T, mainly as the parameters describing the flows of a generating 
system for it. 

In this sense, we introduce the following definition that extends the notion 
of strict nilpotent integrability and encompasses the notion of completely 
integrable systems too. 

Definition 14. A vector field T is said to be complete nilpotent integrable 
if there exists an Abelian algebra f)r of symmetries ofT, with T E f)r, such 
that 

fir 2 N' • 

Let us assume, consistently with our previous assumptions, that T is 
completely nilpotent integrable and that fjp is finitely generated; we denote 
by Yfj, a family of generators. 

Theorem 2. Let T be a complete nilpotent integrable dynamical system with 
a finitely generated Abelian algebra of symmetries fir such that A f C hr - 
Then there exists a family of higher-order constants of the motion gy such 
that 

n 

and conversely. 

Proof. Since T E fir, then T = a^Y^ with E F. Now, for any deriva¬ 
tion Y in the symmetry algebra of T, we have that [Y, T] = 0 implies that 




NILPOTENT INTEGRABILITY AND REDUCTION 


19 


Y (a^) = 0. Hence, the elements are annihilated by all derivations in the 
symmetry algebra of F, in particular by all elements in fir, more specifically 
by all derivations tangent to the projection Y — N ; then they are in TV". □ 


If 9p are local functions parametrizing the flows of X^, then we have that 

r (9J = 9tt , 

and r (g^) are lower order constants of the motion whose evolution will be 
described by Eqs. (19). 

Clearly, this picture reduces exactly to the case of completely integrable 
hamiltonian systems in the case when 

hr 2 C' D M'. 


4. Generalized algebraic reduction of nilpotent integrable 

SYSTEMS 

In this section, we will discuss a generalized reduction procedure, adapted 
to the framework of nilpotent integrability, that extends the well-known re¬ 
duction techniques for vector fields on manifolds. The discussion of general¬ 
ized algebraic reduction will be done as before on the realm of differentiable 
algebras. Later on, particular instances of our construction and the relation 
with well-known reduction theorems will be discussed. 

If we are given a dynamics T on a differentiable manifold M, we can 
reduce it in two ways: first by restricting it to an invariant submanifold and 
secondly, by quotienting it with respect to a suitable projection map. Both 
procedures are perfectly illustrated in the classical example of a hamiltonian 
flow with symmetry, first restricted to an invariant submanifold defined by 
level sets of its constants of the motion and then quotiented with respect to 
the projection map defined by the existence of cyclic variables. 

Hereafter we shall describe these two approaches to the reduction in the 
setting of differentiable algebras, because it makes the approaches indepen¬ 
dent on the existence of auxiliary geometrical structures. 

4.1. Restriction to invariant subspaces. Given a differentiable algebra 
T with spectrum X = Spec K (J-"), an affine differentiable subspace Y C 
X is determined by a closed ideal J of the algebra T. Thus if T is the 
differentiable algebra J~(M ) of smooth functions on the manifold M and 
S C M is a regular submanifold, then the set Js of functions vanishing at S 
defines a closed ideal of -F(M) and F(M)/Js = X{S). In turn, any closed 
ideal J of J- defines a differentiable algebra, given by T / J . Elements in 
the quotient space T / J will be denoted by / + J or simply by [/], unless 
a specific realization in terms of functions on a given set is provided. 

Definition 15. Given a dynamics T on T and a closed ideal J of T, we 
shall say that J is T-invariant ifT(f) E J for any f E J . 
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In this case, the derivation T induces a derivation [r] on the algebra F / J 
by defining [r](/ + J) = T(/) + J . Such a derivation will be said to be 
the restriction of T to F / J . It is clear that if the T-invariant ideal J is 
the ideal of functions vanishing on the submanifold S, then [r] defines a 
dynamics on the submanifold S which is nothing but the restriction of the 
vector field Y to S. 

4.2. Quotienting by invariant relations. The algebraic analog of an 
equivalence relation on a manifold M is a subalgebra of the algebra of func¬ 
tions F. In fact, let us assume that ~ is an equivalence relation on M such 
that the quotient space M = M/ ~ is a manifold and the canonical projec¬ 
tion 7 r: M — 5 - M is a submersion. Then the pull-back with respect to 7r of 
the differentiable algebra of smooth functions on M determines a subalgebra 
7 t*F{M) of F(M). 

Conversely, if F C F ( M ) is a subalgebra, we can define an equivalence 
relation on points of M by declaring that x r^j y iff f(x) = f(y) for all 
/ £ F. Given a differentiable algebra F, not every subalgebra is a differ¬ 
entiable algebra itself, however this is true for closed subalgebras. Then if 
F is a closed subalgebra of F and X = Spec K (T'), X = Spec K (T") are the 
corresponding spectra, there is a natural continuous projection 7r: X —>■ X 
such that 7T* = i is the canonical inclusion of F into F. As in the situation 
above of restriction by an invariant subspace, if we are given a dynamics T 
and a closed subalgebra F : we shall say that it is T-invariant if T(/) 6 F 
for all / £ F. In this case, the restriction of Y to F defines a derivation of 
the subalgebra F. This restriction, denoted by T, will be called the quotient 
of T with respect to the subalgebra F. 

If F = 7 t*F(M), as in the example above, then the vector field T is n- 
projectable and induces a vector field T on M iff F is T-invariant. The 
vector field T is just the quotient of T with respect to ir*F(M). 

Notice that a given equivalence relation, determined by a subalgebra F, 
does not have to be compatible with a given subspace determined by an 
ideal J . We shall say that the equivalence relation F is compatible with 
the subspace J if J C F. 

4.3. Generalized reduction. Both procedures discussed above, restriction 
to a subspace and quotienting by an equivalence relation, can be combined, 
i.e. we can restrict the dynamics to a given subspace and quotienting it out 
by an equivalence relation. In fact we can proceed in two different ways: 
first we restrict with respect to the ideal J and then quotient with respect 
to the closed subalgebra F, or viceversa first we quotient with respect to 
the subalgebra F and then we restrict with respect to a given ideal J . The 
first procedure leads to a dynamics [T] on the algebra F / {F n J) and the 
second gives a dynamics [T] on (F + J)/J. Indeed, the smallest equiva¬ 
lence relation containing F and compatible with the subspace defined by J 
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is the subalgebra T + J. The ideal J is an ideal of F + J, hence T + J / J 
is an algebra. The later reduction in a geometrical form can be traced 
back to the Marsden-Weinstein reduction theorem in the realm of Hamil¬ 
tonian systems with symmetry (Ma74| : the former, in an algebraic form, 
traces back to Grabowski-Marmo reduction theorem for Poisson manifolds 
[Ga94j . Consequently, we will call the first reduced dynamics [T] the geo¬ 
metric Marsden-Weinstein reduction and the second one [T], the algebraic 
Grabowski-Marmo reduction. Both reduction schemes lead to the same 
dynamics as is proven in the following theorem. 

Theorem 3. Let J- be a differentiable algebra with a unit element and T 
be a derivation. Let J be a closed T -invariant ideal and J- be a closed 
T -invariant differentiable subalgebra of T. 

The reduced subalgebras T / J H J7 and JF + J / J are isomorphic differen¬ 
tiable algebras. Moreover, the reduced dynamics are transformed into each 
other by the previous isomorphism, that is [f 1 ] = [T]. 

Proof. The isomorphism among the algebras J n TjT and T + J / J is 
provided by the second isomorphy theorem. We introduce the map 

Hf) = f + J, 

whose kernel is J PI the first isomorphy theorem provides the conclusion 
we were looking for. We will denote the isomorphism so constructed by [<f>]. 

Finally, let us observe that 

[$]*([f])(/ + J) = f (/ + JnJ r ) = f (/) + j n t = [f](/ + j n f). 

□ 

In what follows, since the algebraic Marsden-Weinstein and Grabowski- 
Marmo reduced algebras are isomorphic as well as the corresponding reduced 
dynamics, any of the isomorphic spaces and dynamics obtained in this way 
will be denoted by THjF and the induced dynamics by T j p. 

4.4. Reduction of nilpotent integrable systems. We shall discuss now 
what happens to the integrability properties of a system under algebraic 
reduction, namely how do the constants of the motion and the nilpotent al¬ 
gebras behave under algebraic reduction. We shall focus only on the simpler 
situations, namely the extension to nilpotent integrable systems of Jacobi’s 
procedure of elimination of nodes, Poincare’s reduction of order, or reduction 
by cyclic variables in classical mechanics fSinTH] . 

Consider a T-invariant family Q C Af. If g £ Q, then T(g) E Q , i.e., all 
descendants of the higher-order constant of the motion g must belong to 
Q. Thus the tower T(g ) starting at g must be contained in Q. Hence a 
T-invariant family Q in the nilpotent algebra is the union of towers T(g) of 
higher-order constants of the motion. Because of this, we will also say that 
the family of higher-order constants of the motion Q is nilpotent complete. 
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Let J be a closed ideal generated by a T-invariant family of higher-order 
constants of the motion Q = {g a }. It is clear that such an ideal is T-invariant 
because the tower T(g a ) C Q, i.e. T(g Q ) G Q for any g a G Q. The restricted 
algebra defined by J is given by X / J . 

Let us consider a subalgebra X C X defined by a family of Abelian sym¬ 
metries s C hr, be., / G J 7 if X^(f) = 0 for all G s. We shall denote the 
previous subalgebra as X = X s . 

The following statement follows easily from the theory developed above. 

Theorem 4. Let T be a nilpotent integrable system on the differentiable 
algebra X with Cartan algebra fir- Let Jg = XQ be an ideal generated by 
a nilpotent complete family of higher-order constants of the motion Q, and 
s C fir a family of Abelian symmetries ofT. The reduced dynamical system 
Tjjr defined on the reduced algebra XHj g X s defined by the ideal Jg and the 
subalgebra X s is nilpotent integrable. 


Proof. Due to the fact that we are in the generalized reduction scheme de¬ 
scribed in previous Sections 4. » 


we may use Theorem [3| 


We check first that the reduction of T to the invariant subspace defined by 
the ideal Jg is nilpotent integrable. It is easy to ascertain that the nilpotent 
algebra Mj g of the restriction [T] of T with respect to the ideal Jg contains 
the nilpotent algebra A f restricted to X j Jg. Then 


( 20 ) 


M^ g c(Af/(Mnj g y. 


On the other hand, maximal Abelian subalgebras h [f] of the restricted 
derivation [T] are larger than the restricted Abelian subalgebra [fjp]. Since 
the system is nilpotent integrable, we have M' C fir- Then, due to Eq. (20), 
we get 

(Ab/(Abn Jg)' c [f)r] • 


Hence we conclude that 

M' Jg c [f)r] cf) [r] , 

and the restricted system [T] is nilpotent integrable. 

Consider now the quotient system T defined on the subalgebra F = X s . 
We have that the nilpotent algebra M of T is just A f H X s , i.e., T = T a = 
M n X. At the same time, because T is nilpotent integrable, we deduce that 
A f C l)r, hence 

M' =M' + X' =M'+s". 

Observe that since X = X s , then X is the invariant set s'\ therefore X' = s”. 
But because of the stability Lemma [3j s being Abelian implies that 5 " is 
Abelian too and s C s" . Now, hr is a maximal Abelian algebra of T, so 
5 C hr implies that s" C hr- We conclude that 


Af' c hr = h F , 


and the system is nilpotent integrable. 


□ 
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5. Some explicit classes of nilpotent systems. The 

HIGHER-ORDER CalOGERO-MOSER SYSTEMS 


5.1. Linear systems. A simple instance of derivations is provided by deriva¬ 
tions on C°°(M n ) which are linear with respect to the linear structure on M n . 

If r is a derivation on C°°(R n ), it is linear if [T, A] = 0, where A is the 
derivation whose flow is given by = ehx. In global linear coordinates 

x l on R n , the derivation A is the dilation vector field A = x l d/dx' 1 . If 
T is linear, there exists a linear map A: R n —>■ R n such that T = Xa = 
A J i x l d/dx 3 , where Xa is the derivation whose flow is given by ift{ x ) = 
e tA (x). The algebra C°°(R n ) is generated by any collection of linear functions 
f a (x ) = (a, x), a € (R n )*, containing a generating set oq, .. ., a n for (R n )*. 

A simple computation shows that £p f a = fA*a, where A* is the adjoint 
map to A. Hence the algebra of constants of the motion C for T is gen¬ 
erated by the linear functions f a where a G ker A*. One can see that the 
real spectrum of C is given by ker A and the natural projection is given by 
7 r: R n —>- ker A, defined as ir(x)a = (a,x), for all a € ker A* = (ker A)*. 
Alternatively, we can choose a basis on ,a r of ker A*. Then the map 7r 
can also be written as 7r = f ai v l where v l is the dual basis to a t . 

In a similar way we obtain the nilpotent algebra AT of T as the subalgebra 
generated by the linear functions f a where ( A*) k a = ( A k )*a = 0 for some 
k. Thus if A is nilpotent, then A f = C°°(M n ) and the system is strictly 
nilpotent integrable (and it possesses a generating set of linear higher-order 
constants of the motion). 

In the case of an arbitrary linear system T = Xa, we may use the Jordan 
decomposition of A in its semisimple and nilpotent part, i.e., A = S + N 
where S is diagonalizable, N is nilpotent and [S, A] = 0. Then, the algebra 
generated by the powers of S defines a Cartan subalgebra for T and the 
system is trivially nilpotent integrable. 

5.2. Reduction of higher-order free systems. 


5.2.1. Reduction of the free system of order two: The Calogero-Moser sys¬ 
tem. As in the Introduction, let r (t) = (x(t),y(t), z(t)) be the coordinates 
of a particle moving freely in R 3 , i.e., satisfying r = 0. Then we may define 
a 2 x 2 real symmetric matrix 


*m = ( x ; z 



that satisfies the equivalent system 


( 21 ) 


X = 0. 


As we already discussed, we can reduce this system in order to obtain new 
dynamical systems exhibiting non-trivial potential terms. We can formulate 
abstractly this situation by considering the tangent space M = TTL to the 
space of 2 x 2 real symmetric matrices TL. The points in M will be denoted 
by (X,V), X. V £ TL . Consider now the second order system on TL defined 
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by X = 0. Its flow is given by X(t) = Xq + Vo t, Xq. Vo E %. Notice that 
the vector field corresponding to such system is 

P2) t-v±. 

with the obvious matrix multiplication product understood. 

We may write the matrix X E T~L as: 

(23) X(t) = xl + ycri + za 3 , 

in terms of the matrices k = 1,2,3^} 

'0 A _ / 0 1' 

,1 oj ’ a2 - 1-1 0 , 

The group 17(1) = 50(2), realized by matrices 


CTi = 


O' 3 = 


1 0 
0 -1 


G = 


cos f sin tp 
— sin f cos f , 

acts by conjugation on Ti, X <—> GXG -1 , and transforms any matrix X of 


the form (23) into a diagonal matrix: 
(24) 


X = GQG~ l = G (^ ) G~ l 


We can also write the previous relation as 


X = GM W G 


-l 


with M (0) = Q , 


for later convenience. Let q = <72 — qi denote the difference between the 
eigenvalues of Q. The infinitesimal generator r of 17(1) determines a basis 
for its Lie algebra: 

T = G~ l G = fa 2 . 

Below we list some useful commutation relations: 


(25) [cxi, cr 2 ] = —203 , [02,03] = -2ai , 

(26) [< 7 i, Q\ = qa 2 , [cr 2 , Q] = q<j\ , [a 2 , Q\ = qa 1 , 

(27) [< 72 ,G]= 0 , [ 01 ? G\ = 2 < 7 3 sin p , [ 0 - 3 , G\ = 2ai sin f . 

(28) [r, Q] = qipai, [r, Q] = qpai , [r, ai] = 2pcr 3 . 

Computing the first derivative of a curve X (t) using the factorizacion Eq. 


(24), leads to: 


X = G(Q + [r, <3])G _1 = GM^G- 1 . 

with 

(29) M (1) = Q + [ t ,Q] = M (0) + [r, Af (0) ] . 

Then, we get immediately: 

(30) X = G(M (1) + [t,mW])G _1 = GM^G- 1 


Notice that 02 here is i-times the second Pauli matrix. 





NILPOTENT INTEGRABILITY AND REDUCTION 


25 


with 

(31) M (2) = M (1) + [r, Af (1) ] . 

Iterating the computation, we get 

(32) X (fc+1) = GM^+^G” 1 , Af (fc+1) = M (fc) + [r, M (fc) ], k > 0 . 


Notice that Eq. (30), exhibits explicitly the U( 1) symmetry of the system 
X = 0 and provides its equivalent formulation 

MW = -[r,AfW], 

which is also a Lax representation for it. Equivalently, substituting the 
expression for AfC) given by Eq. (29) we get 

Q = -2qipa 3 + (</3 + 2 qip)ai . 


This relation, by working the commutators out (using Eqs. (26)) and taking 
into account that the matrix Q is diagonal, leads to the system 

(Q= ~2qip 2 a 3 
|0 = q(p + 2qip. 

These three equations can be explicitly written as (using q = q\ + 52 ): 

1= 0 

q = 4 qy 2 


ip = 


2 g . 

— ip. 

q 


The third equation can be easily reduced by a quadrature to 

C 

= ~2 > 
qz 

where C is a constant. In fact, this equation can be viewed as the construc¬ 
tion of an invariant C := q 2, p, and the previous system becomes: 

1= 0 

4C 2 

q = 


(33) 


r 
2 q . 

(p = - ip = 

q 


2 Cq 


However, thinking in the higher-order situation, we prefer to proceed em¬ 
phasizing the role of consta nts of the motion and the generalized reduction 

Clearly the quantity L = [ X , X], represent- 


4.4 


scheme discussed in Section 
ing the angular momentum, is a constant of the motion: 

d 
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In coordinates (q, p) we get: 


(34) 


L = [X, X] = -q pa 2 = -la 2 , l = -tr(Lcr 2 ), 


where t is the third component of the angular momentum, and 


(35) 


1 


l 


La 2 = q pl 2 , p = —-^tr (La 2 ) = 


2 q 2 


T 


We may use any level set of the function l to define an invariant subspace. 
The restriction to it of the vector held T, Eq. (33) reads 

h 2£ 2 

Q = -5-^3 

q 3 

2 iq 

ip =- 3 - • 

<r 

It shows that the vector held d/dp dehnes an Abelian symmetry algebra of 
T. Quotienting with respect to it we finally get, in terms of the coordinates 
qi,q 2 , the classical Calogero-Moser system: 

2l 2 
qi = 


(36) 


<?2 = 


(<?1 - ^ 2 ) 3 ’ 
2£ 2 

(q 2 -qi) 3 ' 


5.2.2. Reduction of the system X = 0: The third-order Calogero-Moser sys¬ 
tem. We will discuss now various reductions of the third order system X = 0 
that represents an uniformly accelerated system. 

By analogy with the previous discussion, the system is dehned in the 
space M = T 2 TL , whose points will be denoted by (X,V,A), X. V, A £ TL. 
The vector held associated with it is 


r 


~ v ~k +A 


d_ 

dV' 


We may certainly perform a reduction of the system by means of the constant 
of the motion X = A = const. However, we prefer to perform a different 
reduction that would mimic the angular reduction discussed before for the 
system X = 0. Besides, in doing so we will show how to proceed to treat 
systematically higher order ‘free’ systems of the form X ^ = 0, k > 3. 
Notice that the system T is trivially nilpotent integrable with how X(t) = 


^0 + 

Keeping with the same notation as before, the recursion equations (32) 
provide 


X = GM (3) G _1 with Af (3) = M (2) + [r, Af (2) ] . 


Thus, the [/(l)-invariance of the system leads to the reduced system 
(37) A/( 2 ) = -[ t ,mW] , 
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with M given by Eq. ( |31[ ). Before computing the explicit expression of 
the system (37), we will discuss the role of higher order constants of the 
motion associated to the angular momentum in this context. 

The angular momentum L = [ X , X] is not anymore a first order constant 
of the motion: 

L = ±[X,X] = [X,X\^0. 


However, 


dt, 


(38) 

and, finally, 

(39) 


L = ^[ X ,X] = [X,X] 


L=^ s {X,X] = 0, 


thus L is a third-order constant of the motion. The quantity L defines 
the tower of higher-order constants of the motion L^ = L, L ^ = L and 
L 1 ' 1 - = L, where the upper index indicates the order of the constant of 
the motion as in Section |3.4[ More generally, if we consider the equation 


X = 0, the quantity L will be a constant of the motion of order 2k — 3 
(except for n = 1 where L = 0). 

As we know from the previous discussion, Eq. (34) gives L = —q 2 ipa 2 - 
Then it is easy to ascertain that all the derivatives Lx* are proportional to 
(72. Then, multiplying by <72 we can consider the quantity which is a 

multiple of the identity. Taking their traces we obtain the scalar functions 


(40) 

h = \ty(Lg 2 ) 

= q 2 v, 

(41) 

£2 = £\ = ^Tr (La 2 ) 

= 2 qqip + q 2 ip , 

(42) 

£3 = £2= *Tr(L<7 2 ) 

= 2 q 2 y + 2 qqip + AqqCp + q 2 


As in the second-order case, £3 is the third component of the generalized 
‘angular momentum = L and £3 = 0. We can integrate the relations 
£2 = £3 ? and £\ = £2 to obtain: 


£\ — 2 ^ 30 ^ ^ 10 ’ 
£2 = £3ot + ^20 5 
£3 = £30 , 


with £10, £20 and £30, are initial values for the quantities £\, £2 and £3 re¬ 
spectively. At this stage, we can use the tower of angular constants of the 
motion £^ to compute M^ k \ First observe that 

M I ' V} = Q + [ t ,Q] = Q + qtpai = Q + —ay, 

q 
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and 

Af (2) = Af (1) + [r, Af (1) ] = Q + 2%a z + —o\ . 

r q 

It can be written, for purposes that will become clear in the subsequent 
computations, in the form 

Af® = Q + Mf V 3 + MfVi, 

with Af® = 2£ 2 /q 5 and AfV = t^/q. Hence we get for Af® the expression 
(43) Af® = Af <( 2 ) + [r, Af®] 

= Q + (a f ® + 20 Af®) as + (0g + Af® - 20 Af ®) <n 
= Q + AffVs + AffVi, 


with 


Af,® = 


Af® = 


Af ® + 20Af ® = 6 


<7 3 


-4 

<r 


m + m< 2) - 20Mf = h + - 4 


<7 


( r 


r 


Hence, as Q is diagonal, the relation Af® = 0 provides the system of equa¬ 
tions 


(44) 


>S r hhq~hq 
Q = -6-4-CT 3 , 

0 = k + Mi-pi - A. 

q q- q b 


Thus we have obtained from the generalized angular momentum reduction 
of the third-order free system, a new integrable system. 

Definition 16. The system of equations in M 2 

Jih{qi - <72) - h(qi ~ <72) 


(45) 


<7i = 6 


q 2 = 6 


(<72 - <7i) 4 

M 2 V 2 - <7i) - h{(l2 ~ <7i) 


(<72 - <7i) 4 

will be called the third-order Calogero-Moser system. 


This system is nilpotent integrable because it has been obtained by re¬ 
duction from the nilpotent integrable system X = 0. 

Other reductions of the original system can be obtained by using a differ¬ 
ent tower of constants of the motion. For instance, we may have used the 
‘energy function’ 

E = -X 2 , 

2 

hence, 

E = X ■ X , E = X 2 , E = 0 . 
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Let 


ei = TV (E) = -(q\ + q 2 ) + q 2 f> . 


Then, using the tower of functions e\, e 2 = e\, e 3 = (e 3 = 0), we 
can obtain different reductions of the original system, all of them nilpotent 
integrable. 

There are further reductions of the third-order Calogero-Moser system 
that finally will take us back to the second-order Calogero-Moser system 
and, in addition, become Hamiltonian. The simplest way to do that is to 
use the constant of the motion A = X. Denoting by c\ = |(du + A22), 
C 2 = A 12 , and C 3 = \{A\\ — A 22 ), we deduce that the Hamiltonian function 

„ 1 / -2 -2i 2V? 


Cl 1 

-y (91+ 92 ) + 2 (92 - 9 i)(-C 2 sin 2 tp + c 3 cos 2</?) 


is a constant of the motion. This Hamiltonian has been obtained from the 
obvious Hamiltonian for the system X = A. namely 

H = ^(x 2 + y 2 + z 2 ) - (cix + c 2 y + c 3 z ), 

using the coordinates qi : q 2 ,<P as before and the angular constants of the 
motion. 

However we cannot reduce this Hamiltonian further, because it is not 
invariant with respect to the Abelian symmetry d/dip, except when C 2 = 
C 3 = 0. In this case, the equations of motion Eq. (45) reduce to the simple 
form 

2i\ 


(46) 


91 = 


92 = 


(91 - 92) 3 
2 i\ 


+ Cl, 


+ Cl. 


(91 - 92P 

which is just the Calogero-Moser system with a constant acceleration term. 


5.2.3. The fourth-order Calogero-Moser system. We will conclude the dis¬ 
cussion on the reduction of higher-order free systems by considering the 
angular reduction of the fourth-order system X = 0. Proceeding again as 
in the third-order case, we get: 

X (4) = GM ( 4 ) G _1 , with M( 4) = M 3 + [r, M (3) ]. 

By using the structure of given by equation Eq. ( |44[ ), we obtain 

M (4) = + M ( 3 3 ) a 3 + M[ 3 ) a 1 + [T,Q] + M^ ) [T,a 3 ] + M[ 3 ) [T,a 1 ] 

= + Mj 4) CT3 + MfV 

with the recursive equations: 

M= M^ ] + 2 (pM[ 3) , 

M[ A) = <p'q+ M{ 4) - 2 ipM^ ] 
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Thus a simple computation taking into account the explicit expressions for 

A/f( 3 ) 


3 and ' given in Eqs. (44) 


M (4) = Q (4) + ^ 


leads us to 
U\ + 8^1 £3 U\q + 324 l 2 q 

q 3 q A 


+ 


8 l\q 2 8£f 


Q 


q‘ 


+ 


U J\ Q - i?Jl 0 hqq~hq 2 1/1^2 lo ?'jq 

-r Z-,-O-o-14—j-io—rr- 

q q z q J q° q° 


This gives us the fourth-order system 
(47) 


q (4 ) = 6^2 + 8£i4 + 4^g + 32£ig 2 g 


8^fg 2 8£ 4 


0"3 


o-i 


6. Discussion and Future Perspectives 

In this work, we have tried to identify the main features which are needed 
to provide an explicit description of the flow of a given dynamical system 
and that do not rely on any additional geometrical structure. The nilpotent 
algebra consisting of the family of all higher-order constants of the motion, 
together with a maximal Abelian algebra of symmetries of the given dy¬ 
namics, are sufficient to realize this description, provided that they satisfy 
a natural relation between them, i.e., that the maximal Abelian algebra 
contains the annihilator of the nilpotent algebra. Natural examples, like 
higher-order free motion, that do not exhibit any obvious geometrical struc¬ 
ture associated to them, fall in this category. We have named these systems 
the nilpotent integrable ones, since their flows can be written as a polynomial 
in the time parameter. 

Such systems generalize in a natural way the notion of Hamiltonian com¬ 
pletely integrable systems. They exhibit canonical forms and the theory of 
reduction applies nicely to them. Also, we have shown how to construct 
examples of them without recurring to any additional geometrical structure 
like a symplectic or a Poisson one. The algebraic language offered by the 
theory of differentiable algebras has proved to be particularly well adapted 
to make manifest the conceptual content of the notion of nilpotent inte- 
grability. It is important to point out that our approach allows to treat 
dynamics on spaces exhibiting mild singularities. 

To avoid to make this paper inordinately large, we do not have discussed 
how the general ideas exposed here particularize when we consider additional 
geometrical structures, or when such structures can be determined from the 
integrability properties of a given system. Neither we have discussed the set 
theoretical counterpart of the theory, that is, we have not made any attempt 
to translate the results contained in this paper in the geometrical language 
of smooth manifolds, vector fields, etc. 

Other issues of primary interest that can be treated in the present theo¬ 
retical framework are for instance superintegrable systems and, in a different 
vein, the extension of the notion of integrability to non-commutative spaces. 
These subjects will all be discussed in subsequent papers. 
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